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Abstract 

We consider scalar-tensor theory for describing varying speed of 
light in a spatially flat FRW space-time. We find some exact solutions 
in the metric and Palatini formalisms. Also we examine the dynamics 
of this theory by dynamical system method assuming a ACDM back- 
ground and we find some exact solutions by considering the character 
of critical points of the theory in both formalisms. We show that for 
any attractor the form of non-minimal coupling coefficient is quadratic 
in terms of the scalar field ^. Also we show that only attractors of the 
de Sitter era satisfy the horizon criteria. 



1 Introduction 

There are some ideas that suggest the constants of nature, such as gravi- 
tational constant and the speed of light, should be time-space dependent. 
[U El E]- Although the number and the role of fundamental constants are 
still debated pE], there are different mechanisms leading to varying constants. 
Brans-Dicke Scalar-tensor theory [21 [5] is the first formulation of a dynamical 
gravitational constant as while as a theoretical explanation of Mach's princi- 
ple is inherent in it. Another important physical constant which has attracted 
considerable attention recently pQ is the speed of light. Theories with varying 
speed of light (VSL) have been firstly proposed by Moffat, Albrecht, Magueijo 
and Barrow [El [7] as an alternative to the inflation mechanism solving some 
problems of Big-Bang cosmological models [61 [8]. In their formulation the 
Lorentz invariance is broken and there is a preferred frame, called cosmolog- 
ical frame, in which the speed of light is only a time dependent field. In this 



1 



frame there exists a pre-set function [7J E] representing the speed of light and 
enters in FRW equations as an input. 

It is a well-known fact that it is possible to have a varying speed of light 
theory and preserving the general covariance and local Lorentz invariance 
|10j . The price that have to be paid for this, is to introduce a time-like 
coordinate x° which is not necessarily equal to ct. In terms of x° and x, one 
has local Lorentz invariance and general covariance. The physical time t, can 
only be defined when dx°/c is integrable. 

The most general scalar-tensor action of gravity which allows for a dy- 
namical speed of light is illustrated in [TT] . This action is previously analyzed 
by many authors using the metric approach. Demianski et al [12] present a 
class of cosmo logical models derived from Noether symmetry requirement. 
These models describe accelerating evolution of a FRW universe filled with 
dust matter and exhibit power-law dependence of coupling and potential to 
the scalar field. There is also some tracking solutions of this model [T3j , in 
which the time evolution of the scalar field tracks the expansion rate of the 
universe. 

In our earlier paper [H] we have found the exact classical cosmological 
solutions assuming an exponential coupling between a scalar field, represent- 
ing dynamical speed of light, and the geometry, with or without cosmological 
constant. Here we shall continue our previous work [H] on investigating the 
exact cosmological solutions with varying speed of light. In the following sec- 
tions we shall use it and find it's exact cosmological solutions for the spatially 
flat universe. We get the results firstly for a metric theory. Then we discuss 
the exact solutions using the Palatini approach in which the connection and 
metric are independent degrees of freedom. In the last section we examine 
the dynamics of this theory by dynamical system method assuming a ACDM 
background. By considering the character of critical points of the theory we 
find some exact cosmological solutions in both formalisms. 

2 The Model 

The Jordan-VSL action which we use here is the one presented in [TT] : 

S = J d A x^g-{F{^)R-2U{^)-Z{^) g^d^d^)+S m [&, g, u ] (1) 

in which F(*f?) = (c/co) 4 and U(^>) are arbitrary regular functions of the 
scalar field representing the coupling of the scalar field \I/ with geometry 
and it's potential energy density respectively. Cq is a constant velocity and 
hereafter we sail put 8nG = c 4 = 1. The first part of the above action func- 
tional is the gravitational part, including Ricci scalar R and a dynamical term 
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for the velocity of light with arbitrary coupling function Z{^>). The latter is 
the action for matter fields, 0j, and doesn't involve the scalar field so the 
matter is minimally coupled to gravity. As emphasized in the introduction, 
here we have assumed that there is a time-like coordinate x°, which is not 
equal to ct and thus dx°/c is not necessarily integrable. The dynamics of ^> 
depends on the functions F and Z. But note that Z can always be set equal 
to unity by a redefinition of the field Therefore only one arbitrary function 
remains. 



2.1 Metric approach 

In the metric approach the metric and the scalar field \I/ are dynamical vari- 
ables. The variation of the action ([T]) with respect to them gives: 

F(V)(R^ - \g^R) = + d^d u m - ^u(d a ^) 2 

+ V^F(^) - g^V^WF^f) - g^U^i) (2) 
ldF dU 

Also the weak equivalence principle holds because the matter fields are min- 
imally coupled to the metric. This implies: 

= (4) 

where the energy-momentum tensor of matter is defined as: 

TP= 2 ^ 
V ~9 ognu 

In a cosmological context, applying the field equations (J2])-(|4j) to FRW uni- 
verse in which the metric has the following form: 

ds 2 = - dx ° 2 + a ( t )2( + r 2 df] 2 } (6) 

leads to the following cosmological equations: 

3F(H 2 + 4) = p + l^ 2 - 3HF + U (7) 
a z 2 



2F{H-^) = (p + p) + ^ 2 + F-HF 

or 
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( ™*) = 3f (* + 2H » + *)-§, (9) 

p + 3# (p + p) =0 (10) 

These are c-variable FRW equations, the wave equation of ^ field and the 
conservation law respectively. H(x°) = \-^o is the Hubble parameter, p and 
p are the energy and pressure densities of a perfect fluid considered as matter 
field and dot denotes derivative with respect to the time-like coordinate x°. 
These equations form a coupled set of nonlinear differential equations for 
H(x°) and *(a?°). The t fine-like coordinate x° is related to cosmic time by 
the relation: 

dt = — =pm ( n ) 

In cosmological application dx° / c is integrable and gives the physical time. 
Therefore the physical Hubble parameter H p (t) = can be evaluated as 
H p = H(x°)^. 

Substituting ^oy^o by in (Tini) gives: 

f t +3H p (p + p) = (12) 

This shows that in this model the conservation equation ( TTTJT) is valid even in 
terms of cosmic time. However it is possible to have the non- conservation of 
energy-momentum if we change our model. The way that the conservation 
relation changes, highly depends on the model. For example in the preferred 
frame approach [6J, violation of the energy-momentum conservation occurs 
because of appearing a source term proportional to the gradient of c in the 
conservation equation. 

The conditions which one should impose on VSL models are usually in- 
spired by the cosmological puzzles. In order to solve the horizon problem of 
the standard cosmology, one should set [lj a/ a — c/c > for the early uni- 
verse and also one has d > . These lead to some constraints on the range 
of possible values of integration constants, appeared in the solutions. 

The cosmological solutions of the greatest interest are those for which 
the time evolution of the Hubble parameter is proportional to the inverse 
of the cosmic time (corresponding to power-law expansion) or a constant 
(corresponding to de Sitter expansion). Considering a spatially flat FRW 
universe, we shall distinguish two cases. A c-dominated universe by which 
we mean S m = U = 0. A (c-A)-dominated universe means that S m = but 
A is not zero and is of gravitational type. 
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2.1.1 c-dominated universe: 

Putting S m = U = in the equations of motion, we get two independent 
equations: 

3FH 2 = ^ty 2 - 3HF (13) 

- 2FH = \i> 2 + F — HF (14) 

Assuming a power-law dependence for the coupling coefficient F, the above 
equations have the following solutions: 

H ~b (15) 

~ x° a (16) 

And also the coupling function F(ty) ~ \l/ 2 which is a particular case emerged 
by requiring the existence of Noether symmetry [12]. The cosmic time is 
defined as (assuming a^2): 

t ~ (x°y- a/2 (17) 

Thus the physical Hubble parameter and the speed of light can be written 
as: 

H p ~j^ a~t u (18) 

c~t^ (19) 

Requesting an expanding universe together with horizon criteria lead to the 
following constraints on the range of possible values of constants: 

v > v> (20) 

2 — a 

Taking now a = 2, the corresponding solutions are given by substituting 
a = 2 in the relations (TIoT) and ([TBI . But now the time is defined as: 

t~lnx° (21) 

so that: 



H p ~ cons. a ~ e vt (22) 



c ~ e K * (23) 

So a = 2 leads to de Sitter expansion for cosmic scale factor. In this case it 
is needed that: 

v > v -k>0 (24) 
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Moreover another solution which leads to the power-law expansion is: 

H~cons. i)~ e ax ° F ~ ip 2 (25) 

This is a special choice which is used previously by many authors [TQ1 HH [T5] . 
The cosmic time is: 

t ~ e~^ x ° (26) 



so: 

H p ~j^ a~t u (27) 
Considering both criteria pointed before, one gets the following constraint: 

v > (28) 

2.1.2 (c — A)-dominated universe: 

A mentioned before, this era corresponds to a matter free universe for which 
the potential U is nonzero and has the form U = AF in which A is a constant. 
Demanding power-law expansion for cosmic scale factor one can easily show 
that the solution is: 

H ~ cons. ip ~ e ax ° F ~ ^ (29) 

And the cosmic time is: 

t ~ e'^ ax ° (30) 
Here a > 0, so the speed of light is decreasing in cosmic time. 

2.2 Palatini approach 

In the Palatini approach the metric and connections are considered as in- 
dependent fields. The variation of the action (CQ) with respect to the metric 
gives: 

F(*)(R^ - l -g^R) = + dfld,,* - l -g^{dJ>) 2 - g, u Um (31) 

where R^ v is Ricci tensor constructed from connections and R = g^ u R^ v . By 
varying the action ([!]) with respect to the connection, one arrives at another 
field equation: 

V«(V=£<rF(*)) = o (32) 
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in which V represents the covariant derivative with respect to the connection. 
This equation shows that the connections are Levi-Civita connections of a 
metric related to as: 



or equivalently: 



V = F{V)g^ (33) 

(34) 




where: 

1 
2 



X| 7 = -(5p 7 ]nF(V) + %d p hiF{*) - g M g a5 d 5 \nF^)) (35) 



is the difference between the affine connection and the Christoffel symboles 



and finally the variation with respect to \I> gives: 



It is worth noting that in this case the matter energy-momentum tensor 
is divergence free with respect to covariant derivative defined with Levi- 
Civita connection of the metric, i.e. V M T^ = 0. This emplies that the test 
particle shall move on the metric geodesic calculated using the Levi-Civita 
connection. An explicit proof of this point for a more general action can be 
found in [16] . 

Equation ( 1331) shows that the F-field or equivalently the c-field acts as a 
conformal factor of space-time metric, g^. Using the transformation rules of 
the Riemann and Ricci tensors under rescaling f[3"3"j) and then inserting those 
in ( !3"T1) and ( l36l) . we obtain: 

i o JP 2 i 

+ V„$,F(tt) - ^V a V a F(^) - ^t/(*) (37) 

VaV °* = ^ " " | V « Vai? (*) + 2^^ F ) 2 ) ( 38 ) 

where = g^R^fJ and R^J is the Ricci tensor constructed from Christoffel 
symboles. In the spatially flat FRW space-time one can easily shown that the 
field equations (j3"7|) and f[3"8j) together the conservation law lead to: 

k vi/ . 3F 2 

mH 2 + - 2 )=p + --3HF- — + U (39) 
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-2F(H+—) = {p + p) + V 2 -HF + F- 



3F 2 



(40) 



2F 




dU 



(41) 



(42) 



As we did for the metric approach, one can find some exact solutions for 
the spatially flat universe in two cases: c-dominated and(c — A)-dominated. 
The interesting point which can be easily shown is that demanding power- 
law expansion for cosmic scale factor in the Palatini framework leads to the 
solutions which are the same as the results found using the metric formalism. 

3 Dynamics of VSL theories 

Another way to find out some exact solutions for cosmological models is the 
dynamical system method [17J. In this method by choosing some appropriate 
variables, one can convert the field equations of the desired theory to a set of 
autonomous differential equations. Then the critical points of the autonomous 
system describe interesting exact solutions. Also this method helps us to 
check the stability of the solutions. Here we consider a class of VSL theories 
which are described by action (JTJ). It is worth noting that this theory is the 
same as to a general scalar tensor theory, but we should take into account 
that the volume element is defined as dx°d 3 x which is different form the 
canonical volume element. So as we have explained before, since in the field 
equations H is not the physical Hubble parameter and derivatives are with 
respect to x° coordinate, the dynamics of VSL theories should be different 
from the scalar tensor theories. Dynamics of a general scalar tensor theory 
in the Jordan frame using metric approach has been considered in [TS] by 
demanding a background ACDM spatially flat cosmology. 

In section 2, to find the exact solutions, we chose F(x°) such that the 
corresponding solution for cosmic scale factor was physically interested. How- 
ever, here we impose a general form for physical Hubble parameter which is 
related to ACDM cosmology, that is: 



where f2 0r = ~ 10 4 , f2 0m = ^ ~ 0.3 and f2 A = 1 — O 0m — f2 0r , then we 

per per 

seek for the corresponding solutions for scalar field. We also assume U(^f) ~ 
F(\l/) n where n is a constant [19]. In order to express equations (1311 . (1401) and 



H p {zf = H 2 [n 0m (l + zf + flor(l + ^ + ^a] 



(43) 
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(|41|) as a dynamical system of the first order differential equations, we first 
write them in dimensionless form as: 

o \l> /2 U F' 2 F' 
1 = Q + Pr + — + — (44) 

m 3FH 2 QF 3FH 2 F 2 F K ' 

IF _Ap j _ V 2 F' H'F' F" 3F' 2 

H FH 2 3FH 2 F F HF F F 2 1 } 

<S>" H'm' ldF r 3H' 9F' 3H'F' 3F" 3F /2 n 

+ ——— + 6 + — + — — + 



F HF F (Nf H 2F 2HF 2F AF 2 

1 dU 3tf' 



FH 2 F 



(46) 



where ' = = ~£p = Wj€jl&- Now we use a set of dimensionless 
phase-space variables x\...x± similar to those introduced in [IB], that is: 

F' U V 2 Pr 

X2 -3FE 2 ' X3 ~6F> X4 -3FH~ 2 (4?) 

Now defining fl m = 3 p^ 2 we write equations (jHJ) and (J4"5"]l as: 

Q m = 1 - x 4 - x 3 - x 2 - xi + ^x\ (48) 

x[ = 3 + 2— + x 4 + 3x 3 - 3x 2 - (2 + —)x x + -x\ (49) 
H HA 

Differentiating x^ and x 2 with respect to In a gives: 

x ' i = Xi [ Xl -4-2 I L\ (50) 
ti 

x' 2 = x 2 [x 1 (l - n) - 2^-] (51) 

where n = Finally differentiating X3 with respect to In a and using (H6l) . 
we have: 

4 = -2(— + 3)x 3 + + x 4 + 5x 3 + (2n - 3)x 2 - 1) - -x? (52) 

H 2 o 

On the other hand one can easily verify that = jf- + \x\, so by substituting 
this relation to equations (|4"9"|) . (I5TJ1) . (loTj) and (1321 we have: 

3 H' H' 
4 = 3 - (- + tt)xi - 3x 2 + 3x 3 + 2-? + x 4 (53) 



4 = x 2 M^ - n) - 2^] (54) 

H 1 1 1 

x' 3 = -2(-? + 3)x 3 + -x x { Xl + x 4 + 4x 3 + (2n - 3)x 2 - 1) - -x? (55) 

1 H' 

< = x 4 [-x 1 -4-2-/] (56) 

Z H p 

These equations describe the cosmological dynamics of the VSL theory in 
Palatini formalism. By a completely similar procedure one can write the 
field equations of the metric formalism as follows: 

3 H' 3 H' 

x[ = -3 - xx(- + -=r) - 3x 2 + ^4 + 2x 3 + -A + 2 TT ( 57 ) 
2 iip 4 iip 

4 = x a [x 1 (i-n)-2^] (58) 

i .ff' jy' i 

x ' 3 = xi [-x 3 ~ 2 ~ Jf + nx ^\ - 6x 3 - 2 ^ Ex 3 - (59) 

1 i?' 

< = x 4 [-x 1 -4-2-/] (60) 

It is important to note that jf- is not always constant so the dynamical 
equations of the metric and Palatini formalisms are not autonomous and we 
can not find the critical points in any regime. On the other hand we know 

H' 

that -jf- for ACDM background is approximately constant in the matter, 
radiation and de Sitter eras and so we can use the dynamical system method 
for these eras [T5] . The critical points and their corresponding eigenvalues are 
shown in Table (00) for the metric formalism and in Table ([2]) for the Palatini 
formalism. 

Considering Table(j21), we see that the only physically interested coordi- 
nates are R5, Mi, M 3 and A 1; for other points X3 < for any n. For 
n < |, A 1; M 3 and R5 are attractors, but demanding positive values for 
X3 imposes another range to n, say —0.93 < n < —0.13. Note that A = 
V281 + 128n 3 + 924n - 732n 2 and C = y/3 + 12n. It should be noted that 
Ri, Mi and Ai correspond to general relativity with constant speed of light. 
Also it is worth noting that these results are different from those obtained in 
[TS] for similar scalar tensor theory. Considering the critical points of Table3 
we see that M' z is an attractor if — | < n < |. Also in the de Sitter era 
A^ for n < I and A' 3 for n > 5 are attractors. The noteworthy feature of 
the Palatini VSL theory is that there exists no attractor critical point in the 
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Table 1: The critical points of the system (I57|) . (1581) . (159]) and (IdT?|) 



Era Cp X2, X3, X4) Eigenvalues 





Ri(Q, 0, 0, 1) 


(1, -1/2, -2, 4) 




i2a(2, 0, -1, 0) 


(1, 2, 1/2, 5-2n) 


Radiation 


it: 3 (4/3, 0, -1/3, 0) 


(5/3, -1/2, 2/3, f - f ) 


= _2 

Hp 


R 4 (-2, 0, -1/3, 0) 


("l,^" 11 , ^ +11 ,3+2n) 




p / 8 (4n 2 -24n+35) 4(2n 2 -9n-2) n x 
K ^2n-1' 3(2n-l) 2 ' 3(2n-l) 2 ' U / 


/ 4 2n+3 10n-29-A 10n-29+A\ 
V2n-l'2n-l' 4-8n ' 4-8n / 



M^O, 0, 0, 0) (3, -1, 43 _^t3) 

Matter M 2 (2, 0, -1, 0) (0, 0, 1,4 - 2n) 

Sc. = -3 M f 6 2(4+n 2 -4n) 2 n 2 -8n-l n\ ( i 3 (2-n)(-3+C) (2-n)(3+C) \ 

Hp 2 3 V2n-l' (2n-l) 2 ' (2n-l) 2 \ '2n-l' 2n-l ' 2n-l / 



Ai(0, 1, 0, 0) 
A 2 (6, 0, -7, 0) 
De Sitter A 3 (2, 0, -1, 0) 

^ = A 4 (4, 0, -3, 0) 



_g -9+V33+96n -9-V33+96n 

(-1,^,^,3 -6n) 
(-3/2, -3, -2, 1 - 2n) 

(3/2, -2, -1, 2-4n ) 



radiation dominated universe and unlike the metric formalism, for all critical 
points, £3 can be a positive quantity. However, in the Palatini formalism we 
have three critical points R' 3 , M[ and A' 3 , which are meaningless because for 
them F'(^f) 7^ while \J/' = 0. It is also important to mention that again 
R' 2 , M' 2 and A[ correspond to general relativity with constant speed of light. 

It is interesting to note that for any solutions which are extracted from 
the analysis of the critical points, the non-minimal coupling coefficient takes 
the form F(^f) ~ (\& — ^o) 2 - F° r an Y regime and critical point, the speed 
of light takes the form c(t) ~ a{t)~ Xl / 4 in both formalisms. The scalar field 
takes the form *&(t) = \t^ Xl ^ + \l/o and the horizon criteria is x\ > 4 in 
the radiation dominated era. In the matter dominated era, the scalar field is 
\I/(t) = \t~ Xl / 3 + \l/o an d the horizon criteria takes the form x\ > 2. On the 
other hand considering the above tables, we see that the critical points in 
the matter and radiation eras do not satisfy the horizon criteria. However, 
in the de Sitter era the horizon criteria is satisfied by any physical critical 
point and the solution for \l> is *ff(t) = \e~ Xlt ^ + ^q. 
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Table 2: The critical points of the system ([53"]) . (1q"4"1) . (155|) and (15"5"]) 



Era Cp (xi, X2, X3, X4) Eigenvalues 





p/ / 8 4n 2 -24n+35 4(2n 2 -9n+10) n \ 
lV2n-l' 3(2n-l) 2 ' 3(2n-l) 2 ' / 


/n 2n+3 lOn-29-A 10n-29+A\ 
\ ' 2n— 1' 4— 8n ' 4— 8n / 




R' 2 (0, 0, 0, 1) 


(-2, 4, 1+ f~\ l 'f^) 
(0, 1/2, 2, 5 - 2n) 


Radiation 


E! 3 (2, 0,0, 0) 


-"p _ _2 


# 4 (4/3, 0, 1/9, 0) 


(0, -1/2, 5/3, f(l-n)) 




# 5 (-2, 0, 2/3, 0) 


(O,^ 11 , ^ 4 +11 ,3 + 2n) 




M[(2, 0, 0, 0) 


(0, 1, -1, 4 - 2n) 


Matter 


M^(0, 0, 0, 0) 


(3, -1, 4*-^) 


= _3 

//p 2 


irU 6 2(4+n 2 -4n) 2(4+n 2 -4n) n N 
iW 3^2n-l' (2n-l) 2 ' (2n-l) 2 ' U J 


/ 3 (2-n)(-3+C) (2-n)(3+CK 
v >2n-l' 2n-l ' 2n-l ' 




A'i(0, 1, 0, 0) 


/ ,1 o -9+V33+96n -9-V33+96n\ 
I ^> 4 ) 4 / 




A' 2 (6, 0, 2, 0) 


(-4, 9 +f^, 9 -f^,3 Qn) 


De Sitter 


A' 3 (2, 0, 0, 0) 


(-3/2, -4, -2, 1 - 2n) 


rr/ 

f p 


A 4 (4, 0, 1, 0) 


(3/2, -4, -1, 2 - 4n ) 



4 Conclusions 

In this letter, we have considered a class of VSL theories which are described 
by the action (TjQ), in both metric and Palatini formalisms. In section 2 we 
have shown that this model exhibits a power-law or de Sitter expansion for 
the cosmic scale factor. In any case applying both formalisms, we have shown 
that F(^f) is a quadratic function in terms of the dynamical scalar field 
This is a particular case emerged by requiring the existence of the Noether 
symmetry for any cosmological point-like lagrangian of a general scalar tensor 
theory |12j . On the other hand, one aspect of our model is if the scalar field 
can be interpreted as dark energy as well as variable speed of light. Many 
authors studied for example the change of fine structure constant based on 
quintessence [20] • To give a look at this point here, note that in the metric 
formalism, using equations (JTj) and (jHJ), we define u e ff as follows: 

P eff & + 2F + 4HF-2U 

U e ff = — = ■ ■ (61) 

n Peff ^ 2 -6HF + 2U 
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The choice of u e ff = 



— 1 leads to the following equation: 



^ 2 - HF + F = (62) 

In the matter-radiation free universe, by using equations (JTJ) , (JSJ) and (16"2"1) we 
have: 

QFHl + 5H F + F - 2U = (63) 

where H = Ho is a constant. By using these equations one can verify that 
the coupling coefficient takes the form Ffif) ~ \I/ 2 in both c-dominated and 
A-dominated universe. As an example of this point, consider a A-dominated 
universe. The equation (I6"3"j) leads to the following solution for the coupling 
function: 

F ~ e ^(5H + v /H 2 +8A)xO (g 4 ) 

And by substituting this result in (1521 we obtain: 



$ „ e =£(5H + y/H*+8A)x° ^ 

So F(^>) ~ \I/ 2 . Similarly, in the Palatini formalism, we define ui e ff as follows: 



p eff & + 2F + 4HF-^-2U 

Ueff = — — = : : ^ (66) 

Peff ^-QHF-^ + 2U 



2F 

Again, u e ff = — 1 leads to the following equation: 

3F 2 

^ 2 -HF + F- — = (67) 
2b 

By considering the modified Friedman equations of this formalism, equations 
( l39j) and (flOl) . in the matter-radiation free universe, one can see that the 
coupling function satisfying the condition (I6"3l . So by using the equations 
(I6"3"j) and (I6"T|) . we see that -F(\I/) ~ ^ 2 in both c-dominated and A-dominated 
universe. So the consistency with Nother symmetry and having constant 
LU e ff, specially u e ff = —1, are both satisfied for F(^) ~ \1/ 2 . But it should 
be stressed here that unlike the usuall ACDM or quintessence models, the 
expression u e ff = — 1 dose not mean necessarily an accelerated universe (i.e. 
Cf > 0). To clarify this point let us combine (CO) and (IE]) in the metric 
formalism (or (1391) and (flOl) in the Palatini formalism) as follows: 

l = -\[0- + 3uJ )P + (! + 3 ^//)Pe//] (68) 
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Converting the derivatives in terms of x° to derivatives with respect to the 
cosmic time, for a matter-radiation free universe, we have: 



Note that this equation is correct in both formalisms with only different 
expression of p e ff so the following results are common for them. From this 
equation we see that the evolution of the speed of light has a crucial role 
in constructing an inflationary universe. If the speed of light is decreasing 
in cosmic time, then the necessary condition for > is tu e ff < — ~, but 
this is not sufficient condition. On the other hand, if the speed of light is 
increasing in time, then the sufficient condition for Cf > is u e ff < — ~, but 
this is not necessary condition. In the other word, by a constant uJ e ff, the 
sign of Cf can change with time. For a simple example, consider the solution 
(1641) . where H(x°) = H and cj e // = —1- By using equation (1691) we have: 



where a is a positive constant. 

We recall that the action (CQ) which is used in this paper, allows for a 
dynamical gravitational constant as well as dynamical speed of light. The 
only difference is that we should take x° coordinate as a time-like coordinate 
in the latter case. So the equations of variable speed of light and variable 
gravitational constant scalar-tensor cosmology are not similar when we write 
them in terms of the physical time and physical Hubble parameter. To clarify 
better, starting action (CQ) as a variable gravitational constant theory and then 
applying it to cosmology, considering a matter-radiation free universe, it can 
be easily shown that the first term in equation (1691) is absent in this case. 
In this theory uj e ff < — | leads to an inflationary universe if the matter- 
radiation density is negligible. 

In section 3 we have considered the dynamics of the VSL theories in 
order to find out other exact solutions. We have constructed the cosmological 
dynamical system which is constrained to obey the ACDM cosmic history. 
Also by considering the corresponding critical points, we have shown that for 

x l 

both formalisms, variable speed of light takes the form c ~ a~~ in each era. 
But the time-dependence of dynamical scalar field is different in the radiation 

all x-^ 

and matter eras and it has the forms \I/ ~ t~ * and \1/ ~ t s~, respectively. 

By calculating the variables Xi, ...,x 4 for exact solutions of the c- domi- 
nated universe which have been obtained in section 2 and comparing them 
with the results of the section 3, one can show that they are not critical 
points in metric formalism. But by choosing appropriate values for a,v and 




(69) 




(70) 
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k in the relations ffT8"]) - f[2"T)) . they can correspond to the critical points R'±, 
R' b , A' 2 and in Palatini formalism. 

It is worth to note that for any solution which are extracted from the 
analysis of critical points, the non-minimal coupling coefficient takes the form 
F(^) ~ — \l/o) 2 - Also unlike the metric approach, there is no attractor 
critical point in the radiation dominated era for the Palatini approach. It is 
also important to note that critical points in the matter and radiation eras 
do not satisfy the horizon criteria. However in the de Sitter era the horizon 
criteria is satisfied by any physical critical point and the scalar field has the 
exponential dependence to cosmic time. 
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